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Abstract 

In  this  paper  an  analysis  of  remote  atmospheric  magnetometry  concepts,  using  molecular 
oxygen  as  the  paramagnetic  species  is  considered.  The  objective  is  to  use  these  mechanisms  for 
the  remote  detection  of  underwater  and  underground  objects.  We  use  the  coupled  Maxwell-density 
matrix  system  of  equations  to  describe  the  magnetization  of  the  paramagnetic  species  in  the 
presence  of  an  intense  modulated  laser  pulse  and  ambient  magnetic  field.  The  magnetic  dipole 
transition  line  that  is  eonsidered  is  the  -  X^'Z~  transition  band  of  O2  near  762  nm.  The  pump 

is  taken  to  be  a  high-intensity,  pulsed,  polarized  titanium-doped  sapphire  laser.  The  model  is  used 
to  investigate  a  magnetic  anomaly  detection  signature  using  the  polarization  rotation  of  the 
Wakefield  behind  the  pump  pulse.  The  major  challenges  are  the  collisional  dephasing  of  the 
atmospheric  oxygen  transitions  and  the  strength  of  the  effective  magnetic  dipole  interaction. 
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I.  Introduction 

Optical  magnetometry  is  a  highly  sensitive  method  for  measuring  small  variations  in 
magnetic  fields  [1-3].  The  development  of  a  remote  optical  magnetometry  system  would  have 
important  applications  for  the  detection  of  underwater  and  underground  objects  that  perturb  the 
local  ambient  magnetic  field.  In  our  remote  atmospheric  optical  magnetometry  model,  a  high- 
intensity  laser  pulse  is  employed  to  drive  the  magnetic  spin  precession  in  the  presence  of  the 
earth’s  magnetic  field.  Zeeman  splitting  of  the  molecular  energy  levels  and  consequently  the 
precessional  frequency  are  proportional  to  the  local  magnetic  field.  This  can,  in  principle, 
provide  a  precise  means  to  measure  variations  in  the  Earth’s  magnetic  field.  For  a  number  of 
magnetic  anomaly  detection  (MAD)  applications,  //G  magnetic  field  variations  must  be  detected 
at  standoff  distances  of  approximately  one  kilometer  from  the  sensor  [4]. 

In  this  paper  we  consider  atmospheric  molecular  oxygen  as  the  paramagnetic  species  in  a 
nonlinear  remote  optical  magnetometry  configuration  depicted  Fig.  1 .  The  propagation  of  the 
high-intensity  pump  laser  pulse  to  remote  (~  km)  detection  sites  is  considered.  By  employing  the 
optical  Kerr  effect  we  show  that  high  laser  intensities  (10^^  W/cm^)  can  be  propagated  to  remote 
locations.  Using  a  high-intensity  polarized  laser  pulse  we  consider  the  polarization  rotation  of 
the  Wakefield  behind  the  pump  pulse.  The  pump  intensity  is  limited  (<  10^^  W/cm^)  to  avoid 
photoionization  processes. 


High  intensity, 
polarized 
laser  pulse 


Water/Ground 


Figure  1:  Remote  nonlinear  optical  magnetometry  configuration.  The  Earth’s  magnetic  field  is 
«  0.5  G  and  5B  is  the  perturbation  caused  by  the  underwater/underground  object. 


Molecular  oxygen’s  paramagnetic  response  is  due  to  two  unpaired  valance  electrons. 
The  ground  state  of  oxygen  ,  commonly  referred  to  as  “triplet  oxygen,”  has  total  angular 
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momentum  /  =  1 ,  total  spin  5  =  1  and  three  degenerate  sublevels.  The  excited  upper  state  being 
considered  is  denoted  by  .  It  has  7=0  and  is  a  spin  singlet  5=0  state  with  only  one 

sublevel.  The  upper  state  can  undergo  three  radiative  transitions,  “  -1)  ’ 

=  0) ,  but  the  latter  is  insignificant  because  it  is  an  electric  quadruple  transition. 
There  is  an  intermediate  state,  referred  to  as  ,  into  which  the  excited  O,  molecule  can 

<5  ^ 

decay  and  is  discussed  in  Appendix  A.  The  O2  transition  line  being  considered  is  the 

-  X^'L~  transition  band  of  oxygen  near  762  nm.  In  the  low  laser  intensity,  long  pulse, 

regime  this  transition  has  been  investigated  theoretically  [5,6]  and  experimentally  [7]  and  is  a 
prominent  feature  of  air  glow. 

In  this  paper,  a  high  intensity,  polarized  titanium-doped  sapphire  laser  is  considered  for 
the  pump  laser.  These  lasers  have  an  extremely  large  tuning  range  from  660  nm  to  1180  nm,  and 
can  have  linewidths  that  are  transform  limited.  A  major  challenge  for  this,  as  well  as  any  remote 
atmospheric  optical  magnetometry  concept,  is  collisional  dephasing  (elastic  collisions)  of  the 
transitions.  The  elastic  molecular  collision  frequency,  at  standard  temperature  and  pressure 
(STP),  is  =3.5x10^  sec”^ ,  where  cr  is  the  molecular  cross  section  and  is  the 

thermal  velocity  [7].  On  the  other  hand,  the  Larmor  frequency  (spin  precession)  in  the  Earth’s 
magnetic  field,  »  0.5G ,  is  O  ^  =  qB^  /  (2mc)«  4.5 x  10® rad / sec  ( TiQ.^  =3 x  10“^eV),  where 
m  and  q  are  the  electron  mass  and  charge  and  c  is  the  speed  of  light.  Since  the  dephasing 
frequency  is  far  greater  than  the  Larmor  frequency,  the  parameters  are  somewhat  restrictive  for 
remote  atmospheric  magnetometry.  However,  rotational  magnetometry  experiments  based  on 
molecular  oxygen  at  STP  and  magnetic  fields  of  ~  10  G  have  shown  measurable  linear  Faraday 
rotational  effects  [7]. 

II.  Atmospheric  Propagation  of  Intense  Laser  Pulses  (Focusing  &  Compression) 

The  magnetic  anomaly  detection  concepts  considered  here  rely  on  propagation  of  intense 
laser  beams  in  the  atmosphere.  Atmospheric  propagation  of  intense,  short  pulse  lasers  are 
strongly  affected  by  various  interrelated  linear  and  nonlinear  processes  [8].  These  include 
diffraction,  Kerr  focusing,  group  velocity  dispersion,  spectral  broadening  and  self-phase 
modulation.  Pulse  compression  can  be  achieved  by  introducing  a  frequency  chirp  on  the  pulse; 
however,  for  the  problem  under  consideration,  pulse  compression  is  not  significant.  Nonlinear 
transverse  focusing  results  from  the  Kerr  effect.  In  general,  a  laser  pulse  propagating  in  air  can 
be  longitudinally  and  transversely  focused  at  remote  distances  (>  km)  to  reach  high  intensities 
(~10^^  W/cm^),  as  indicated  in  Fig.  2. 
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Figure  2:  Simultaneous  transverse  focusing  and  longitudinal  compression  of  a  chirped  laser  pulse 
in  air  due  to  nonlinear  self-focusing  and  group  velocity  dispersion. 

Here  we  present  the  model  describing  longitudinal  and  transverse  compression  of  a 
chirped  laser  pulses  in  air  [8].  For  optimally  chosen  parameters,  the  longitudinal  and  transverse 
focal  distances  can  be  made  to  coincide  resulting  in  a  rapid  intensity  increase  near  the  focal 
region. 

The  laser  electric  field  is  given  by  E(r,  rj,  r)  =  (1  /  2)E{r,  rj,  -i-  c.c. ,  where  E  is  the 

complex  amplitude,  O)  is  the  frequency,  r  is  the  radial  coordinate,  x  =  t-zl c  and  r]  =  z  are  the 
transformed  coordinates,  the  propagation  distance  z  and  time  t  are  in  the  laboratory  frame. 
Substituting  this  field  representation  into  the  wave  equation  results  in  a  higher  order  paraxial 

wave  equation  for  E(r,7,r)  [8], 


(1) 


where  the  wavenumber  is  k  =  col  c  and  the  group  velocity  dispersion  parameter  is 

=  2.2  X 10“^^ sec^  /  cm .  The  nonlinear  refractive  index  of  air  is  n  =  \  +  rij^I ,  where  the  Kerr 

nonlinear  index  is  =  3xl0“'®cm^  /  W  and  7  is  the  laser  intensity.  The  values  for  and  tif. 

apply  to  air  at  STP  and  ?i  =  2k Ik  ~ 762 nm . 

Equation  (1)  can  be  solved  by  assuming  the  pulse  is  described  by  a  form  that  depends  on 
certain  spatially-dependent  parameters.  With  this  assumption,  a  set  of  simplified  coupled 
equations  can  be  derived  for  the  evolution  of  the  spot  size,  pulse  duration,  amplitude  and  phase 
of  the  laser  field.  Taking  the  laser  pulse  to  have  a  Gaussian  shape  in  both  the  transverse  and 
longitudinal  directions,  the  complex  amplitude  can  be  written  as 


E{r,ri,x)  =  EXny' 


(2) 
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where  E^{ri)  is  the  field  amplitude,  6{r])  is  the  phase,  Rirj)  is  the  spot  size,  a{ri)  is  related  to 
the  curvature  of  the  wavefront,  T{r])  is  the  laser  pulse  duration,  and  Pin)  is  the  chirp  parameter. 
The  quantities  ,  6,  T,R,  a,  P  are  real  functions  of  the  propagation  distance  rj .  The 
instantaneous  frequency  spread  along  the  pulse,  i.e.,  chirp,  is  5co{ri,  r)  =  2p{rf)T  /  T^irf) ,  where 
P{t])  =  T{t])l{2p,)dTit])ldri  <  0(>  O)  results  in  a  negative  (positive)  frequency  chirp,  i.e., 

frequency  decreases  (increases)  towards  the  back  of  the  pulse. 

Substituting  Eq.  (2)  into  Eq.  (1)  and  equating  like  powers  of  r  and  z  ,  the  following 
coupled  equations  for  R  and  T  are  obtained. 


d^R 


drj^  k^R^ 


1- 


PmlTj 


d^T 


4A  £o 


1 


k  P,,R}T^ 


+ 


API 


(3a) 


(3b) 


where  6^  =  P(0)r (0)  is  proportional  to  the  laser  pulse  energy  and  is  independent  of  rj , 

P{rp  =  nRi^  {rplirj) !  2  is  the  laser  power,  7(77)  (7)7(8;^)  =  /(0)/?^(0)r(0)/(/?^(7)r(;7))  is 

the  peak  intensity  and  7^^  -  ^  !  {pnrifP  is  the  self-focusing  power.  In  Eq.  (3)  the  initial 
conditions  are  given  by  a(0)  =-(kR(0)  /  2)dR(0)  /  dr/  and  p(0)  =  T(0)  /  (2ppdT(0)  75^  =  0. 

The  first  term  on  the  right  hand  side  of  Eq.  (3a)  describes  vacuum  diffraction  while  the  second 
term  describes  nonlinear  self-focusing,  i.e.,  due  to  .  Nonlinear  self-focusing  dominates 

diffraction  resulting  in  filamentation  when  P  >  P^^  w3GW  [8,9]. 

In  the  limit  that  the  pulse  length  does  not  change  appreciably  the  laser  spot  size  is  given 
by  R(j/)  =  R(0)  [1  -  2a(0)j/  7  Z^q  -i-  (1  -  P  7  -i-  (0))(77  7  p  ^  ,  where  Z^^  =  kR^  (0)  7  2  is  the 

Rayleigh  length.  The  spot  size  reaches  a  focus  in  a  distance  r/ 7 Z^q  =  a(0) /(l-P/  7^^  -1- (0)^ 


as  long  as  P  <  (1  -1-  (0))7^ 


NL  * 


Eigures  3(a)  and  3(b)  show  the  evolution  of  the  laser  spot  size  and  the  intensity  as  a 
function  of  propagation  distance  for  X  -  762nm .  At  focus,  the  laser  intensity 


(^focus  =6xlO^°W7cm^)  and  spot  size  {R^^^  =1.3  mm)  are  held  constant  by  choosing 
appropriate  initial  conditions:  wavefront  curvature  a(0)  =  37 ,  pulse  duration  r(0)  =  100  psec , 
and  chirp  yff (O)  =  0.  By  changing  the  laser  energy  and  the  initial  spot  size,  the  nonlinear  self- 

focusing  effect  changes  the  focal  point  from  0.25  km  to  0.75  km  (see  Eig.  3).  Nonlinear  laser 
pulse  propagation  allows  for  moving  of  the  detection  site  location. 
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Figure  3:  Evolution  of  (a)  laser  spot  size  and  (b)  normalized  peak  laser  intensity  as  functions  of 
propagation  distance  for  different  initial  laser  energies  and  spot  sizes.  The  laser  energy  and 
initial  spot  size  for  the  solid,  dashed,  and  dotted  lines  are  =  100,150,190  mJ  and 
/?(0)  =  4.7, 6.7, 8.2  cm ,  respectively. 

III.  Optical  Magnetometry  Model 

The  four  levels  of  O2  being  considered  in  the  magnetometry  model  are  shown  in  Fig.  4. 
The  ground  state  is  split  by  the  Zeeman  effect  into  three  levels  1 1) ,  1 2) ,  and  1 3)  and  the  excited 
state  is  denoted  by  1 4) .  The  magnetic  quantum  number  m  associated  with  the  various  levels  is 
indicated  in  Fig.  4.  The  excited  state,  level  1 4) ,  can  be  populated  by  right  hand  polarized  (RHP) 

light  from  level  1 3)  or  by  left  hand  polarized  light  (LHP)  light  from  level  1 1) .  Here  the 

quantization  axis  and  the  direction  of  the  static  magnetic  field  are  taken  to  be  along  the  direction 
of  laser  propagation,  z  -  axis.  Circularly  polarized  radiation  carries  angular  momentum  ±h 
which  is  directed  along  the  propagation  direction,  hence  the  selection  rule  for  allowed  transitions 
is  Am  =  ±1 ,  which  will  conserve  angular  momentum  [10].  It  should  be  noted  that  this  transition 
is  strictly  magnetic  dipole-  and  spin-forbidden,  but  spin-orbit  coupling  between  the  and 

(m  =  0)  states  leads  to  a  transition  with  a  magnetic  dipole-like  nature  and  a  larger  than 
expected  dipole  moment  [5,6,11]. 
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Figure  4:  Energy  levels  associated  with  the  ground  and  excited  state  of  O^.  The  transition 
frequency  corresponds  to  ~  762  nm  =  1.63eV) .  The  Zeeman  splitting  of  the  ground  state  is 
caused  by  the  ambient  magnetic  field. 

A  high-intensity  pump  pulse  generates  a  magnetization  current  density  =  cVxM , 

where  M  is  the  magnetization  field.  The  current  density  in  turn  generates  a  response  electric 
field  and  can  also  modify  the  pump  pulse.  The  response  electric  field  E  is  given  by  (Gaussian 

units)  - (1  / c^)d^  /dt^jE  =  (4/r / c^)dj ^  I dt  =  (Att / c)d(V x M) / dt .  The  magnetization  is 

represented  by  a  sum  of  RHP  and  LHP  components  M(z,r)  =  -i-M^(z,r)e^  -hc.c., 

where  (z,r)  =  N (z,r) ,  =  N p^p^^izj) ,  N  is  the  density  of  the  paramagnetic 

species  (oxygen  molecules),  is  the  effective  magnetic  dipole  moment  associated  with  the 
transitions,  p^^  and  p^^  are  the  off-diagonal  coherence  of  the  allowed  density  matrix  elements 
(see  Fig.  4)  and  =  (e^  +iey)  /  2  are  vectors  denoting  the  polarization  direction.  The 

magnetization  current  density  can  be  written  as  7^  =-/c5M^  (z,r)/5ze^  -I- 
ic dM^  ( z,  r)  /  -I-  c.c.  In  terms  of  the  x  and  y  components. 


=-/(c/2)5(M^-MJ/5ze^  +i(c/2)5(M^ +MJ/5ze^ +c.c. 


The  density  matrix  equation  is  given  by  [10,12,13]  dp^ldt=  -ico^  p^  + 
ni  Pin  1  rclaxation  terms ,  where  denotes  the  interaction 


frequency,  the  phenomenological  relaxation  terms  are  due  to  elastic  and  inelastic  collisions,  and 
spontaneous  transitions  and  the  magnetic  dipole  interaction  Hamiltonian  is  •  B  (Appendix 
B).  The  off-diagonal  coherence  elements  of  the  density  matrix  for  the  relevant  transitions  ( 
1 1)  ^  1 4)  and  1 3)  ^  1 4) )  are  given  by 
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(4a) 

(4b) 


5/741  ldt=- +  iQ.^,  {p,^-p^)  +  /Q43 Ai ^ 

5/743  !  ~ycPA3  “(^43/^43  “*‘(^43  (/^33  ~ /^44  )  "*”  ^^41  A3  ’ 

where  y^  is  the  elastic  collision  frequency  (not  population  transferring),  the  full  set  of  density 
matrix  equations  are  given  in  Appendix  B. 

The  pump  laser  field,  which  induces  the  magnetization  field,  is  expressed  as  a  sum  of 
RHP  and  LHP  fields  Bp^j^p(z,0  =  S^(z,Oe^  +B^(z,Oe^  +c.c. ,  where  ^^^(z,?)  = 
and  \i/{z,t)  =  kz  —  o)t.  The  interaction  frequencies  associated  with  the  allowed  transitions  are 
Q43  =  p^Bj^{z,t) /  fi  =  and  ^41  =  p^B^iyz,t) /  fi  =  Q^(z,0e'’^^^’'^ ,  where 

^RL=  Pm^R  z,  /  ^  is  (half)  the  Rabi  frequency  associated  with  the  RHP  and  LHP  components  of 

the  pump.  Note  that  the  Rabi  frequency  is  defined  with  respect  to  the  peak  field. 

Although  we  are  considering  a  magnetic  dipole  transition  it  is  convenient  to  express  the 
Rabi  frequency  normalized  to  an  electric  dipole  moment.  The  magnitude  of  the  Rabi  frequency 

can  be  written  as  =  iP,n  I A )  (A^peak  >  =  ( A  '  A )  ( A  >  ^)(8^^  /  ,  where 

I  =  cE^  /  (8;r)  is  the  pump  laser  intensity  and  is  the  peak  field.  Taking  the 

normalizing  electric  dipole  moment  to  be  p^=qrR=  2.5x10“^^  statC-cm,  where  is  the 

Bohr  radius  the  magnitude  of  the  Rabi  frequency  is  ^j^^,,j[rad/sec]= 

2.5x10*  {p^  ///^)-^/[W/cm^] .  As  an  example,  for  /=10"W/cm^  and  p^l p^  =10^,  the 

Rabi  frequency  is  =  8x10^  rad /sec. 

IV.  Faraday  Rotation  of  Wakefields  Driven  by  Intense  Laser  Pulses 

The  incident  pump  field  is  taken  to  be  polarized  in  the  x-direction 

E  =  .E^(z,t)A'*^“®'^(e^  +e^)  +  c.c.,  where  o)  is  the  carrier  laser  frequency  and  the  complex 
pulse  amplitude  E^ (z, t)  can  be  modulated.  Employing  the  variables  T  =  t-zl c  and  ri  =  z, 

E  (e^  +6^)  +  c.c. ,  the  corresponding  magnetic  field  in  the  y-direction  is 

B  =  -/E^(r)A'®^c'^''(e^  -e^)+c.c. ,  where  tsk^k-col  2nikN  /  (hy^')  is  the 
wavenumber  mismatch  (damping)  obtained  from  the  linear  dispersion  relation.  The 
characteristic  wavenumber  mismatch  for  X  =  762nm  at  atmospheric  molecular  oxygen  density 
V  =  5.7  X 10'*  cm“*  and  an  equilibrium  population  of  /7jj  =  1  / 3  is  iA^  =  -T^  =-1.7  x  10“*cm'^ 

(1  /  «  6  m) .  As  the  pulse  propagates  through  the  atmosphere  it  induces  a  magnetization 

current  which  generates  a  field  polarized  in  both  the  x  and  y  directions.  The  wave  equation  for 
the  forward  propagating,  y-component  of  the  complex  field  amplitude  is 
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{d ! dri+iNc)Ey{r],T)  =  —iljrN ,  where  the  Faraday  rotated  field  is 
E  {r],  e  +  c.c.  The  magnetization  current  is  a  function  of  the  off-diagonal 

coherence  terms  of  the  density  matrix  elements  p^^(ri,T)^p^.^{T)e~'‘“''  and 

■  The  slowly  varying  quantities  P4^{t)  and  p^^ir)  are  given  by 
reduced  density  matrix  equations  (5/ dr  -ils.(X)^^) p^j^{t)=  and 

{d! ^T-i^co^^)p^^{T)=iil^^{T)p^,v^htvQ  Q^jCr)  =  -«x4(r)/^ ,  =ip^E^{T)lh, 

Po  =  Al  =  A2  =  As  =  1  /  3  ,  A4  =  0’  +  ^Tc  ^  «43  =  «41  =  ^nd  it 

has  been  assumed  that  c |M|  /  o  □  1 . 

In  the  case  where  the  pump  amplitude  is  not  modulated  in  amplitude  or  frequency,  i.e., 
conventional  Faraday  rotation  within  a  long  pump  duration,  5/5r  =  0,  the  spatial  change  in  the 
Faraday  rotated  field  is  given  by  (d!  dr] + iAk)E  (rj,  r)  =  IrikN p^  (E^  /  K) p^Q. ^  ly^ .  After 
propagating  a  distance  L,  the  ratio  between  the  Faraday  rotated  and  incident  intensities  is 
Iy/Io  =  \Eyf/\ES  =  i27r)\LIA)\NplpJhf  {kljylf  . 

In  the  present  model  the  pump  pulse  consists  of  a  pulse  train,  as  shown  in  Fig.  5,  in  which 
the  duration  of  the  individual  micropulses,  denoted  by  ,  can  be  comparable  or  longer  than  the 

damping  time  1  /  .  However,  the  time  separation  between  the  micropulses  T  is  taken  to  be 

long  compared  to  a  damping  time.  With  this  ordering  of  timescales  the  individual  pump  pulses 
excite  the  density  matrix  elements  p^^  and  p^^  which  generate  a  wakefield  that  decays  behind 

the  individual  pump  pulses  (Fig.  5).  The  magnetization  currents  generate  a  Faraday  rotated 
wakefield  behind  each  pump  pulse. 
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Pump  Pulse  Train  & 
Wakefield 


T  »  1/Yc 


Figure  5:  Pump  pulse  train  and  induced  polarization  rotated  wakefields. 

The  general  form  of  the  off-diagonal  coherence  elements  is 

T 

(d/dT-iAa)^Jp^jT)=i^^^  (t)p^  with  solution  p„„(T)  =  ip^^dT  (r')exp(-/Ary„^(r' -r)) 

0 

within  the  pump  pulse.  The  solution  behind  the  pump  pulse  is 

P43  (j)  =  As  ( A  )exp(/Ary43  (r  -  ^  ))  and  p^^  (r)  =  p^^  (r^  )exp(/A(»4j  (r  -  A  ))  .  The  reduced  wave 
equations  for  the  x  and  y  components  of  the  wakefields  are 

(a  /  a/7  +  (n,  r)  =  -(2^  /  c)  J^(T)  =  -  CJcp^^  {t)E^  ,  (5a) 

(a  /  a^  +  iAk)E^  ip,  t)  =  -{In  /  c)  J  (r)  =  i  CJip^^  (r)4 ,  (5b) 

where  k  =  colc,  /y»|a/ar|,  c|A^|  and  Cg—27riNpl^lK)lY^^6'x.\Q~^  isaunitless 
parameter.  In  estimating  we  have  taken  the  magnetic  dipole  moment  to  equal 
Pm  —  Pe  ^10^  =  2.5x10”^^  statC-cm,  the  collision  frequency  to  be  Yc  —  3.5x10®  sec”'  and  the  O2 
density  to  be  W  =  5.7  xlO'^cm'^.  The  current  densities  are  J^xiT)  =  (NpmCo/ 2)ip^^iT)  - p^^it)) 
and  (t)  =  -  iNpmp) /  2)( A3  (r)  +  Aj  (r)) ,  behind  the  pulse  ( r  >  )  they  are  given  by 


WJj). 


h 


n^+Y 

O  /  c 


■.^pIe^Po 


h 


-W(T), 


(6a) 

(6b) 
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where  the  time  dependence  of  the  wakefield  is  captured  by 

W^{r)  =  [cos(Q„(r - cos(Q„r) - 

(Q„  /7j(sin(Q^(r-rp))-e"^'''’  sin(Q^r))  , 

W^{t)  =  [sin(Q^(r  -  r^))  -  sin(Q^r)  + 

(^o !  Yc )  (cos(Q„  (r  -  ))  -  cos(Q„r)) 

1 

0.5 
0 

0  12  3  4  5 

y,  (x-Xp) 

X  10~^ 

1 

^  0.5 
0 

0  1  2  3  4  5 

y  (x-x ) 

c  p 

Figure  6:  (a)  x-component  and  (b)  y-components  of  the  wakefield  response  functions  W  (t) 
behind  the  pulse  for  =  4.5  x  10*’  rad  /  sec .  The  pulse  durations  for  the  solid,  dashed,  and  dotted 
lines  are  =0.1  ns,  0.5  ns,  1  ns,  respectively. 

Figure  6  shows  the  wakefield  time  dependence,  Eqs.  (7),  for  pump  pulse  durations  of 
Tp  =  0.1,  0.5,  and  1  nsec,  pump  pulse  energy  of  100  mJ,  and  spot  size  of  1  mm.  This  results  in  a 

range  of  pump  intensities  from  6x10*^  W/cm^  to  6x10^°  W/  cm^ .  Equations  (5)  indicates  that 
E^^IEq  is  proportional  to  ^  (r) ,  if  M  is  neglected.  For  the  parameters  in  Fig.  6,  the 

/V  /V  /V  /V 

normalized  peak  wakefield  amplitudes  are  E^IE^  w  0.5,  1.5,  and  1.6  and  E^IE^  wlxlO  , 

1 .2  X 10”^ ,  and  2x10“^.  There  is  a  tradeoff  between  driving  the  wakefields  with  a  higher 


(7a) 

(7b) 
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intensity  pump  {E^  ~  )  versus  driving  it  for  a  longer  duration  ( ~  )•  As  a  result,  for 


Jf.3' 


Tp>?)l  the  Wakefield  amplitude  begins  monotonically  decreasing 

For  remote  magnetic  anomaly  detection  small  spatial  differences  in  the  magnetic  field 
must  be  measured.  Here  we  consider  measuring  the  differences  in  wakefield  intensities  at  two 
nearby  locations.  The  locations  are  referred  to  as  (1)  and  (2)  and  have  local  magnetic  fields 

and  B^+5B .  The  intensity  of  the  wakefield  at  location  (1)  and  (2)  is  /j  and  ,  respectively. 
The  fractional  change  in  intensity  of  the  y-polarized  wakefield  is 

|/j  -  /jI  /  /j  =  |(J/|  /  /j  w  2  SEy  /  E^^  ,  where  E^^  is  the  y-component  of  the  wakefield  amplitude 

and  5Ey  is  the  difference  in  the  wakefield  amplitudes  between  the  two  locations.  Figure  7 
shows  the  fractional  wakefield  intensities  for  various  values  of  5B .  For  the  values  shown, 


(2) 


X  10  ^ 


y  (r-T  ) 

V  p 


Figure  7:  Fractional  difference  in  the  wakefield  intensity  for  various  fractional  differences  in  the 
magnetic  field  SB  /  B^.  The  pump  pulse  has  duration  x^  =500  psec,  spot  size  1  mm,  and  energy 

100  mJ.  The  differences  in  magnetic  field  and  corresponding  intensities  are  from  two  nearby 
locations. 


The  pump  pulse  energy  is  I{7x/2)R^Xp,  where  R  is  the  spot  size.  For  a  pulse  of  duration 
Xp  —  SOOpsec  «  2/  ,  R  =  0.1cm  and  intensity  I  =  10'°W/cm^ ,  the  pump  pulse  energy  is 

80mJ/pulse .  For  a  pulse  train,  rep-rated  at  =  IkHz ,  the  average  pump  laser  power  is 
{p)  =  f^i(7x/2)R\=mw. 
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It  is  worth  noting  that  at  sufficiently  high  intensities,  the  upper  level,  level  |4) ,  can  be 
populated  resulting  in  a  laser  induced  florescence  signal  to  lower  energy  levels,  i.e.  levels  1 1) 
and  1 3) .  This  process  is  known  as  the  Hanle  effect  and  is  briefly  discussed  in  Appendix  C.  The 
magnetization  current  resulting  from  the  induced  florescence  of  an  x  polarized  pump  laser  is 
Jm  cos(ft>r)[cos(f2^r)e^  -sin(f2^r)ey  J  [10].  Using  polarization  filters,  the  intensity  on 

a  detector  due  to  the  x  and  y  components  of  the  current  density  can  be  measured  separately. 
Taking  the  ratio  of  the  intensities  from  the  x  and  y  components  of  gives  4//y  oc.  cot^(D.Qx). 
Note  that  the  ratio  is  independent  of  the  collision  rate  as  long  as  the  individual  intensities  are 
greater  than  the  inherent  intensity  fluctuations. 

V.  Discussion  and  Concluding  Remarks 

Remote  magnetometry  has  important  applications,  such  as  detection  of  underwater  and 
underground  objects.  Detection  of  the  magnetic  anomaly  caused  by  such  an  object  is  important 
to  Navy  missions.  In  the  laboratory,  under  controlled  environment,  conventional  magnetometry 
techniques  can  be  used  to  measure  extremely  small  magnetic  field  perturbations  (~  juG ). 
Limitations  on  remote  detection  include  collisional  dephasing  of  magnetic  fluctuations,  air 
turbulence,  etc. 

The  paramagnetic  species  considered  here  is  the  oxygen  molecule  which  has  the 
magnetic  dipole  transition  ( )  near  762  nm.  We  considered  an  intense  pump  laser  to 

induce  a  polarization  rotation  of  the  wakefield.  This  transition  is  resonantly  driven  by  a  linearly- 
polarized  pump  laser  pulse.  Our  examples  suggest  that  the  intensity  of  the  rotated  component  of 
the  wakefield  can  be  measured. 

Numerous  issues  remain  to  be  considered,  these  include  signal-to-noise  ratio  limitations, 
magnetic  field  orientation,  pump  laser  absorption  in  the  atmosphere,  detection  of  the  wakefield 
signal,  possibility  of  inducing  O2  lasing. 

Acknowledgement:  We  would  like  to  acknowledge  Dr.  S.  Potashnik  for  useful  discussions.  This 
work  was  supported  by  the  Naval  Research  Laboratory  Base  Fund. 
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Appendix  A.  Transitions  in  Oxygen  Molecule 

Oxygen’s  abundance  in  the  earth’s  atmosphere,  approximately  2 1  % 

A  =  5.7  X 10^^  cm“^ )  and  its  paramagnetic  response  make  it  a  possible  candidate  species  for  a 

remote  optical  magnetometer  [4-7].  Molecular  oxygen  O2  has  two  unpaired  electrons  in  the 
upper  level  of  the  ground  state,  giving  it  a  paramagnetic  response.  The  ground  state  of  oxygen 
X^'Z~ ,  commonly  referred  to  as  “triplet  oxygen,’’  has  total  spin  5  =  1  and  three  degenerate 
sublevels.  In  atmospheric  conditions  near  the  surface  of  the  earth  (pressure  P  =  l  atm ,  total 
number  density  =  2.7  x  10^®  cm“^ ,  and  temperature  T  —  23.5  meV ),  the  ground  state  is  fully 
populated  because  the  next  excited  electronic  state’s  energy,  =  0.98 eV  is  much  greater  than 
the  thermal  energy. 

The  electronic  configuration  of  molecular  oxygen  is  shown  in  Fig.  8.  The  first  excited 
electronic  state  of  oxygen,  ,  is  referred  to  as  “singlet  oxygen’’  and  only  has  one  spin  state 

(5,m)  =  (0,0) .  This  state  has  an  energy  of  =0.98eV,  can  undergo  spontaneous 
emission  via  a  magnetic  dipole  transition  to  the  ground  state  O2  ((^'A^  -  )  or  a  -  X .  The 

a-X  transition  has  a  wavelength  of  1.27 //m.  This  transition  is  dominantly  due  to  the  orbital 
angular  momentum  L  and  has  spontaneous  emission  rate  of  =2x10“^  sec"'  [14]. 

The  second  excited  state  of  oxygen  will  be  referred  to  as  the  upper  state.  It  is  also  a 
spin  singlet  state  with  only  one  sublevel.  The  upper  state  can  undergo  three  radiative  transitions; 
(m  =  ±1) ,  -  X^Z”  {m  =  O)  and  £»'Z^  - a'A^  where  the  first  and  second 

transitions  are  between  the  different  magnetic  sublevels  of  the  ground  state  and  are  referred  to  as 
the  A  band  [11].  The  transitions  will  be  referred  to  as  b-X,l,  b-X,0 ,  and  b-a  respectively. 

The  b -X  transitions  have  an  energy  of  Ey^_^  =  1.63eV ,  wavelength  =  762nm,  and 
frequency  =  2.5  x  lO'^  rad  •  sec"' .  The  calculated  spontaneous  emission  rates  of  the  b-X,l 
and  b-X,0  transitions  are  =0.087 sec"'  and  A|j_x,o  =  1-6 x  10"’ sec"'  respectively  [14]. 
The  radiation  from  the  b-X,l  transition  can  be  seen  in  air- glow,  night-glow  and  aurorae  [14]. 
The  b-X,l  transition  is  magnetic  dipole-  and  spin-forbidden  and  it  is  dominant  over  the  b-a 
and  b-X,0  transitions  which  are  electric  quadrupole  transitions  [11].  This  can  be  explained  by 
a  large  spin-orbit  coupling  between  the  (?'Z^  state  and  the  X^Z"  (m  =  O)  state.  The  spin-orbit 
coupling  results  in  a  mixing  of  the  levels  and  the  b  -  X,  1 .  The  b-a  transition  has  an  energy  of 
E^^_^  =0.65eV  ,  wavelength  =1.9//m  ,  frequency  =  9.9  x  lO'"*  rad  •  sec"'  and 
spontaneous  emission  rate  of  =  1.4  x  10"^  sec"'  [14].  The  transitions  from  the  upper  &'Z^ 
can  be  monitored  by  luminescence  spectroscopy  [14]. 
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Figure  8:  Electron  occupancy  energy  levels  of  Oj  as  two  oxygen  molecules  are  brought  together. 


Appendix  B.  Density  Matrix  Equations 

Interaction  of  an  oxygen  molecule  with  radiation  is  governed  by  Schrodinger’s  equation 
iWxj/ !  dt  =  Hy/  where  H  =  is  the  full  Hamiltonian,  is  the  electronic 

Hamiltonian  after  Zeeman  splitting,  and  is  the  magnetic  dipole  interaction  energy. 

The  wavefunction  ^(r,t )  =  {t)u^  (r)  can  be  decomposed  into  the  orthogonal  energy 

n 

eigenstates  of  O2 ,  (r) .  The  probability  amplitudes  C„  (t)  are  related  to  the  density  matrix 

elements  (t)  =  Q  (?)C'2  (t) .  The  macroscopic  electromagnetic  fields  are  driven  by  a 

statistical  ensemble  of  molecules,  not  a  single  molecule,  and  therefore  it  is  necessary  to  use  the 
density  matrix  equations  and  to  introduce  phenomenological  relaxations  terms,  i.e., 

nm  ^  Pnm  +  Am  -^/m  Pin }  +  felaxation  terms . 

e 

In  our  model  molecular  oxygen  is  treated  as  a  closed  four  level  atom  composed  of  the 
ground  state  Oj  ( )  and  the  upper  level  O2  )  •  The  ground  state  has  three  spin 

sublevels  m  =  -l,0,+l  which  are  referred  to  as  states  |l),  |2),and  |3)  respectively.  The  excited 
upper  level  is  referred  to  as  state  1 4) .  The  complete  set  of  coupled  equations  for  the  density 
matrix  elements,  assuming  a  closed  system,  are  given  by 

5/3ji  !  dt  =  ~(rj2  '^^2\P22  "*”^31^33  ^41  Am  “^41^14)’  (Bl) 

dp22 1  dt  =  Y  i2Pn  ^  21  ^  22)  P22  ^  32A33  ^  42At4  ’  (B2) 
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5/^33  !  dt  =  r^3/^ll  +  ^23P22  “(^31  +  ^ ^2  )  As  ^43/^44  +  /  (^34/743  “  Q43/734  )  ,  (B3) 

—  “  (r4i  +  r42  +  r43  ^  +  /  (^f24^/7^4  —  0^4/74^  ^  +  /  —  Q34/943  ^  ,  (B4) 

^Pa\  !  —  “/41A1  “  ^*^41  Al  "*”  ^  ^41  (Ai  “  A4  )  "*”  ^*  ^43 Al  ’  (B5) 

^As  /  ~  “As  As  ~  ^^43  As  ^  ^43  ( As  ~  A4  )  ^  ^41  As  ’  (B6) 

^P\i>  ^  y xiPii  ^^yiPyi  ^  ^14/^43  ^  ^43/^14  ’  (BV) 


The  population  level  of  state  |  n)  is  given  by  while  the  coherence  between  the  states  are 
given  by  =  p*^  .  The  transition  frequencies  are  defined  as  =<^m~^n  where  %<x)^  is 
energy  of  the  n*  state.  For  example,  the  state  frequencies  are  co^  =  -Q^ ,  <»2  =  0 ,  rUj  =  ,  and 

CD^  =  co^  and  the  transition  frequencies  are  =a>^+Q.^,  =  -20.^ ,  and  =  co^-Q.^, 

where  is  )  transition  frequency  in  the  absence  of  a  magnetic  field.  The 

Larmor  frequency  is  given  by  Q.^  =  qB^  /  (2wc)  where  q  is  the  electric  charge,  is  the  static 
background  magnetic  field,  and  m  is  the  electron’s  mass.  Equations  (B1)-(B7)  imply 
conservation  of  population  levels,  i.e.,  d{p^^  +  P22  +P33  +  A4)  /  =  0  (closed  system).  The 

populations  are  additionally  normalized  unity,  i.e.  Tr(yC>)  =  1 .  The  interaction  frequency  between 
states  m  and  state  n  is  =^„m-  Specifically,  Q.^^(z,t)  =  p^B,^{z,t)  I  h  and 
^41  (^4)  =  ^  where  is  the  effective  magnetic  dipole  moment  between  triplet 

oxygen  and  the  upper  state  and  5^  ^  corresponds  to  is  the  right,  handed  polarization  of  the  pump 
field.  The  rate  equation  for  p^^  considered  since  it  does  not  couple  to  the  those  in 
Eqs.(Bl)-(B7). 

The  rates  and  consist  of  contributions  from  i)  elastic  collisions  (soft,  dephasing 

collisions  with  no  population  transfers),  ii)  inelastic  collisions  (population  transferring)  and 
spontaneous  emission.  The  elastic  collision  rate  is  taken  to  be  the  dominate  rate  and  we  set 
/41  =^43  =  /31  =7c  •  the  absence  of  the  pump  field  and  at  equilibrium  we  have 

Pn  -  P22  -  A3  -Po  P44  =0-  This  implies  that  r2j  =rj2,  Tj,  =^3  and  r23  =r32  and  we 
take  these  rates,  which  include  inelastic  collisions  and  spontaneous  emission,  to  equal  T^ .  In 
addition,  the  rates  T^^ ,  r42  and  T^j  consist  of  inelastic  collisions  and  spontaneous  emission  and 
we  take  these  rates  to  be  equal  to  Tj, .  Taking  the  inelastic  collision  rates  to  be  equal,  i.e. 

T^  =  Tj^ ,  the  density  matrix  equations  become 
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dp, ,ldt  =  -r,  ( Ai  - Pn )  +  i (^4 Ai  -  ^^41  A4  )  ’  (B8) 

SPii  /  ^  =  -To  ip22  -  P22 )  ’  (B9) 

5a33  ! dt  =  —y gip^^—  A33 )  +  ^ ( ^34 A43  ~ ^43^34 ) ’  (BIO) 

^A(4  ^  ^  ~~Y o^PaA  “Atl)  "*"^'(^41^14  (^43^34  “^34/^43)’  (BH) 

^PaX  =  -YcPaX  -^xPaX  +«^41  (Ai  -P44)  +  ^^43Ai  ’  (B12) 

^AtS  !  —  ~YcPa3  ~^^AzPa^  ■*■^^43  (As  ~  A44)  ^^4iA3’  (B13) 

dp,^  Idt  =  -YgP,^ -io\^Pn  +i^xAPA^ -i^A^PxA  ■  (B14) 


The  phenomenological  inelastic  damping  rate  is  given  by  y^  =  3r^  =  3rj,  « 10®  sec  ^  [6].  The 
equilibrium  populations  for  the  ground  state  are  aT  =  P22  ^  Pll  =1/3  and  for  the  upper  state 

pZ-^- 

Appendix  C.  Resonant  Fluorescent  Excitation  (Hanie  effect) 

At  sufficiently  high  intensities  laser  induced  fluorescence,  i.e.  Hanie  effect,  can  be 
considered.  The  Hanie  effect  refers  to  the  depolarization  of  resonant  fluorescence  lines  by  an 
external  magnetic  field  [1,2,10].  It  provides  a  sensitive  experimental  technique  for  a  number  of 
measurements,  including  remote  measurement  of  planetary  magnetic  fields  [15]  and  spontaneous 
emission  rates  [10],  and  spin  depolarization  rates  [16].  It  is  also  the  basis  of  one  of  the  most 
sensitive  methods  for  measuring  the  lifetime  of  excited  levels  of  atoms  and  molecules  [17].  In 
the  presence  of  a  magnetic  field  the  Zeeman  sublevels  of  the  ground  state  are  split,  resulting  in  a 
difference  in  the  resonance  frequencies  for  LHP  and  RHP  light.  The  resulting  phase  difference 
between  LHP  and  RHP  light  which  is  dependent  on  the  ambient  magnetic  field,  alters  the 
polarization  of  fluorescing  radiation. 

To  discuss  this  mechanism  in  more  detail  we  consider  a  short  intense  laser  pulse 
polarized  in  the  x-direction  =  Eg{z')e''‘^  (e^  +Cl)  +  This  is  just  one  of  many 

orientations  and  configurations  of  the  pump  polarization  and  magnetic  field  direction  in  which 
the  Hanie  effect  can  occur. 

The  pump  pulse  is  intense  enough  to  excite  level  1 4)  at  the  expense  of  levels  1 1)  and  1 3) . 
The  pump  pulse  duration  Tp  is  short  compared  to  the  collision  time  which  in  turn  is  short 
compared  to  a  Larmor  period.  As  the  short  duration,  high-intensity  polarized  pump  pulse  sweeps 
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by  it  leaves  behind  an  excited  state  which  fluoresces  with  polarization  components  different  then 
that  of  the  pump.  The  fluorescence  from  the  excited  state  P44  to  states  pn  and  P33  is  described 
by  the  off-diagonal  coherence  of  the  molecular  density  matrix  elements 
P43  =  -inRoTpp44e“'<^"43-!rc)T  andp4i  =  where 

Qro  =-i/J^E^  Ih,  Qrq  =  ijU^E^  Ifi,  co^^  =  co-Q.^  and  co^^  =  co  +  €l^.  The  magnetization  left 
behind  the  pump  pulse  is  M  =  —M^  e~  (e~  ‘  —  e” '  )  -1-  c.c.  where 

=  N {Eg  /  .  The  associated  current  density  is 

=-M^e”^"*^|^ru43|cos((y43r)-/sin((y43r)|e^  -i-ft>4j  |cos(ft>4jr)-/sin(ft>4jr)|e^  J-i-c.c.,  where 
r043,(Pn  □  Yc  •  The  current  density  has  components  in  the  x  and  y  directions  [10]  which,  for 
coU  Qg,  are  given  by 

=-2Mg  e“’'‘^''ft>cos(ft>r)[cos(f2^r)e^  - sin(Q^r)ey  ] .  (Cl) 

By  using  polarizer  filters,  the  time  average  intensity  on  a  detector  due  to  the  x  and  y  components 
of  the  current  density  can  be  measured  separately.  Taking  the  ratio  of  the  intensities  from  the  x 
and  y  components  of  Jpf  gives /^.//y  oc  cot^(floT).  The  ratio  is  independent  of  the  collision  rate 
as  long  as  the  individual  intensities  are  greater  than  the  inherent  intensity  fluctuations. 


18 


References 

[1]  D.  Budker,  W.  Gawlik,  D.F.  Kimball,  S.M.  Rochwester,  V.V.  Yashchuk  and  A.  Weis,  Rev. 
Mod.  Phys.  74,  1154  (2002). 

[2]  Optical  Magnetometry,  D.  Budker  and  D.F.J.  Kimball  (eds.)  (Cambridge  U.  Press, 
Cambridge,  UK,  2013). 

[3]  G.  Bison,  R.  Wynands,  and  A.  Weis,  Appl.  Phys.  B  Lasers  Opt.  76,  325  (2003). 

[4]  J.P.  Davis,  M.B.  Rankin,  L.C.  Bobb,  C.  Giranda,  M.J.  Squicciarini,  “REMAS  Source  Book,” 
Mission  and  Avionics  Tech.  Dept.,  Naval  Air  Development  Center  (1989). 

[5]  R.J.  Brecha,  L.M.  Pedrotti  and  D.  Krause,  J.  Opt.  Soc.  Am.  B  14,  1921  (1997). 

[6]  W.  Happer  and  A.C.  Tam,  “Remote  Laser  Pumping  of  Molecular  Oxygen  in  the  Atmosphere 
for  Magnetic  Field  Measurements,”  Naval  Air  Development  Center  Report  N62269-77-M-4202 
(1977). 

[7]  R.  J.  Brecha,  Appl.  Opt.  37,  4834  (1998). 

[8]  P.  Sprangle,  J.R.  Penano  and  B.  Hafizi,  Phys.  Rev.  E  66,  046418  1-21  (2002). 

[9]  A.  Couairon  and  A.  Mysyrowicz,  Phys.  Rep.  441,  47  (2007). 

[10]  P.W.  Milonni  and  J.H.  Eberly,  Laser  Physics  (Wiley,  Hoboken,  NJ,  2010). 

[11]  B.  Minaev,  O.  Vahtras,  and  H.  Agren,  Chem.  Phys.  208,  299  (1996). 

[12]  R.W.  Boyd,  Nonlinear  Optics  (Elsevier,  Burlington,  MA,  2008). 

[13]  M.  Scully  and  S.  Zubairy,  Quantum  Optics  (Cambridge  U.  Press,  Cambridge,  UK,  1997). 

[14]  B.  E.  Minaev  and  H.  Agren,  J.  Chem.  Soc.,  Earaday  Transactions,  93,  2231  (1997). 

[15]  C.  Kumar,  L.  Klein,  and  M.  Giraud,  in  Work.  Adv.  Technol.  Planet.  Instruments,  edited  by 
J.  Appleby  (Lunar  and  Planetary  Institute,  Eairfax,  Virginia,  1993),  p.  14. 

[16]  S.  Curry,  W.  Happer,  A.  Tam,  and  T.  Yabuzaki,  Phys.  Rev.  Lett.  40,  67  (1978). 

[17]  A.  Corney,  Atomic  and  Laser  Spectroscopy,  (Oxford  Univ.  Press,  Oxford,  2006)  p.  478. 


19 


